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The mean-field dynamics of a Bose gas is shown to break down at time th = (ci /y) In N where 7 
is the Lyapunov exponent of the mean-field theory, N is the number of bosons, and ci is a system- 
dependent constant. The breakdown time tj, is essentially the Ehrenfest time that characterizes the 
breakdown of the correspondence between classical and quantum dynamics. This breakdown can be 
well described by a quantum fidelity defined for one-particle reduced density matrices. Our results 
are obtained with the formalism in particle-number phase space and are illustrated with a triple-well 
model. The logarithmic quantum-classical correspondence time may be verified experimentally with 
Bose-Einstein condensates. 


I. INTRODUCTION 

The nonlinear Gross-Pitaevskii equation (GPE), as a 
mean-field theory, has been the dominant tool in describ¬ 
ing the dynamics of Bose-Einstein condensates (BECs) 
in ultracold atomic gases mm- However, we face a 
quandary when the mean-field dynamics of a BEC be¬ 
comes dynamically unstable or chaotic mm : on one 
hand, one may regard this instability as an unphysical ar¬ 
tifact resulted from the mean-field approximation, since 
the exact dynamics of a BEC is governed by the many- 
body Schrodinger equation, which is linear and thus does 
not allow chaos; on the other hand, the dynamical insta¬ 
bility was observed in experiments [10IH5) and it has been 
proved with mathematical rigor that the GPE describes 
correctly not only the ground state but also the dynam¬ 
ics of a BEC in the large N limit {N is the number of 
bosons) [T51 [T7] . 

Our aim in this work is to resolve this fundamental 
dilemma. Our study shows that the mean-field theory 
(the GPE) is only valid up to time 

Tft, = — InfV -I- o(ln7V), (1) 

7 

where 7 is the Lyapunov exponent of the mean-field dy¬ 
namics and Cl is a constant that depends only on systems. 
With this time scale, the dilemma is resolved: on one 
hand, in the large N limit {N —>■ 00 ), Th goes to infinity 
and thus the GPE is always valid just as proved rigor¬ 
ously in Ref. [HI; on the other hand, the time r/, increases 
with N only logarithmically and it is not a long time for 
a typical BEC experiment. Eor example, for the system 
studied in Ref. [3], the Lyapunov time = I /7 ~ 1 ms. 
As the number of atoms in a BEC prepared in a typical 
experiment is around 10^, we have Th ~ 10 ms. As a 
result, the dynamical instability or the breakdown of the 
mean-field dynamics can be easily observed in a typical 
experiment as reported in Ref. |12j . 
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This time scale Th is essentially the Ehrenfest time, 
which is the time that the correspondence between the 
classical and quantum dynamics breaks down [iHiiin]. 
The usual Ehrenfest time tek = (ci/ 7 ) ln(A//i), where 
7 is the Lyapunov exponent of the classical motion and 
A is a typical action [13]. The similarity is due to that 
the GPE can be regarded as a classical equation in the 
large N limit [20]. Therefore, our result paves a way to 
experimental investigation of a fundamental relation in 
the quantum-classical correspondence — the logarithmic 
behavior of the Ehrenfest time — as N can be varied in 
experiments. 

We cast the quantum dynamics onto the particle- 
number phase space (PNPS), which is a rearrangement 
of Fock states. In this phase space, for a nearly coher¬ 
ent state and in the large N limit, quantum many-body 
dynamics is equivalent to an ensemble of mean-field dy¬ 
namics. When the mean-field motion is regular, mean- 
field trajectories will stay together and the Bose gas re¬ 
mains coherent. If the mean-field motion is unstable or 
chaotic, mean-field trajectories will separate soon from 
each other exponentially, leading to decoherence of Bose 
gas and breakdown of the mean-field theory. So, there 
are two distinct types of quantum dynamics, whose dif¬ 
ference can be characterized by the quantum fidelity for 
one-particle reduced density matrices. 

We investigate the Ehrenfest breakdown numerically 
in the system of a BEG in a triple-well potential pTH35] , 
which may be the simplest BEC model that embraces 
chaotic mean-field dynamics. With this model, we verify 
numerically the Ehrenfest time and show that our quan¬ 
tum fidelity can well capture the characteristics of two 
different types of quantum dynamics. 

The mean-field instability or breakdown has been dis¬ 
cussed in literature [3 EHZl [MSI]- However, a general 
and explicit relation between mean-field chaos, number 
of particles and breakdown time is still lacking. And in 
PNPS not only such breakdown can be understood intu¬ 
itively and quantitatively, but the significance of a local 
phase structure is also apparent, distortion of which leads 
to decoherence. 
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II. PARTICLE-NUMBER PHASE SPACE 

In Ref. [20j . it is shown that many quantum systems 
become classical in the large N limit. A dilute Bose gas 
belongs to this class of quantum systems: its dynamics 
becomes classical and it is well described by the mean- 
field GPE in the large N limit. In this section, we in¬ 
troduce PNPS, where this quantum-classical correspon¬ 
dence in the large N limit becomes transparent. 

A. Definition 

Any quantum state of a system of N identical 
bosons with M single-particle states can be regarded as 
a wavefunction ip{x) over an {M — I)-dimensional lat¬ 
tice space, which we call particle number phase space 
(PNPS), via 


And the wavefunction ip{x) in PNPS has a phase struc¬ 
ture. For any x and y in PNPS, 

M 

arg (f{x) - arg (p{y) = y^) arg (mod 27r), 

(7) 

which shows a wavevector k: ki = N aigipi oc N. This 
phase structure is important as it will give us an estimate 
of the time in our later discussion. It is worth noting 
that when N ^ oo, there is no limit of the wavefunction 
'^{x) because its wavevector k diverges. 

Overall, we find that the coherent state corresponds to 
a single-peaked wavepacket with plane-wave phase struc¬ 
ture in PNPS. In the following, we shall discuss quantum 
dynamics in PNPS and its relation to the mean-field dy¬ 
namics. Note that the formalism of PNPS was also used 
in other contexts [SH [33] , where phase structure and dy¬ 
namics, however, were not discussed. 


M 


'Nxi 


ip{x) = (0| 




1 ^), 


( 2 ) 


where XiS are entries of the M-dimensional vector a:, 
Nxi e {0, ...,A} for I < * < M and '^^Xi = 1. And 
a\ and on are the creation and annihilation operators for 
the Ath single-particle state, with [ai,aj] = Sij and fii = 
alcti- The continuous limit of PNPS is a hyperplane in 
[0,1]^^ (defined by constraint = 1 )j where we 

can define (for i from 1 to M) 

{Xi) = J dxx^\^p{x)\'^ (3) 

((Ax^f} = J dx (xi - (xi}fl(p(x)l^ (4) 

to characterize the average position and spread of the 
distribution |(/j(a;)p over PNPS, given normalized. Of 
course for any finite JV, the integral should be interpreted 
as summations over all x in PNPS. 

As an example of our particular interest, we examine 
an SU(M) coherent state |4')c in PNPS: 





(5) 


where = 1. In such case, we say \ip) (an M- 

dimensional vector with ijji as its entries) is the mean-field 
state of the SU(M) coherent state |'I')^. It is straightfor¬ 
ward to show for this coherent state |4') 


B. Dynamics 


Consider a quite general Hamiltonian of a Bose gas 

M 


H = ^ a|aj -I- 




where Hfj = Hj* and Uij = Uji. Corresponding to the 
Schrodinger equation idt |'I') = H I'k), there is anequa- 


tion of motion (EOM) for ip{x]t) in PNPS (Eq. (Al) in 


the Appendix). We are especially interested in the dy¬ 
namics of a nearly coherent state (p{x;t), which satisfies 
the following two conditions: 

(z) the distribution \(p{x-,t)\‘^ is localized such that 
1/N \J ((Axi)^) ^ 1 for all z = 1,2,..., M; 

(zz) a local wavevector k{x;t) = Va, arg (/^(a;; t) exists 
in PNPS and varies insignificantly over a scale of 1/N, 
i.e., \dxikj \ N for all z, j = I, 2,..., M. 

With these two conditions and keeping only finite 
terms in the large N limit, an approximate (to 0(1)) 
EOM for ip{x-,t) in PNPS can be derived (see Eq. (A2) 
in the Appendix). Mathematically, there are 5-function 
solutions to this EOM (Eq. (A2)): 


M 


Vj(a;; t) = exp[za(t)] n 5{xi - x/{t)) exp[zfc°(t)xi]. (9) 
2=1 

In these 5-function solutions, x/{t), k/{t) satisfy the fol¬ 
lowing equation 


idtp = [Hmf, p ], 


( 10 ) 


{x.) = \H\ ((Aa:.)') = |zA,|^(l-|z/-,n/IV, (6) 


where HuF^jit) = -I- 2UijPij{t) and 


which indicate that the coherent state corresponds 
to a localized distribution \ip{x)'\^ in PNPS that peaks 
around (|z/’iP, |z/^ 2 p,' ’ ’ : IV’Mp) with a vanishing spread 
at large N. 


p,,{t) = ^x0(t)cr0(t)e*('=-(‘)-'=?(*))/'^ . (11) 

This is just the mean-field EOM for the one-particle re¬ 
duced density matrix. 
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Conditions (i) and {ii) reflect our expectations of 
nearly coherent states (see Eqs. Q and 0). The ex¬ 
istence of (^-function solutions corresponds to the estab¬ 
lished result that for any time to, when N ^ oo, coherent 
states at < = 0 stay coherent when t = to m- 

The results above can be interpreted as follows: at 
large N, for any initial state satisfying the two conditions, 
its time evolution may be regarded as the superposition of 
mean-field dynamics of i5-functions, since any function in 
PNPS can be decomposed into a superposition of a cloud 
of (5-functions! This is similar to the quantum dynamics 
of a single-particle wavepacket in real space: it can be 
regarded as a cloud of classical particles and each of them 
follows the Newton’s EOM. 

As the quantum-classical correspondence between a 
quantum wavepacket and a classical particle will break 
down at the Ehrenfest time, the correspondence between 
one state in PNPS and its mean-field description — one 
(5-function solution (see Eq. 0 ) — will also fail when the 
mean-field trajectories of the (5-functions in the cloud di¬ 
verge. 

The breakdown time Th can be estimated using a con¬ 
ventional strategy in quantum chaos as in Ref. |19j . Es¬ 
sentially, before the breakdown the wavepacket of nearly 
coherent states in PNPS expands in the form of exp jt, 
where 7 is the Lyapunov exponent of the mean-field dy¬ 
namics. According to Eq. 0 , for t < Th, 


1 It 


And there is a consistent mean-field description only if 
local wavevectors across the wavepacket are almost equal, 
that is, 

KtN~^A{t) <^1, (13) 


where kN~^ is the average rate of growth of curvature 
dxikj and the N dependence is written explicitly. Sub¬ 
stituting (12) into (13), we have 


qt-I-Int-I-InK ^ (A-I-^) In Af. (14) 


The Ehrenfest time in Eq. 0 is obtained with 
Cl = A -|- i, which is independent of N or 7 . Numeri¬ 
cal verification of this relation will be presented later. 

Note that it is well-known that the quantum-classical 
correspondence may last far beyond the Ehrenfest time 
(see, e.g.. Ref. [31] )• Similarly, it is possible that the 
mean-field theory remains valid even after our first es¬ 
timate t = Th', this interesting and special topic will be 
left for future study. 


III. EXAMPLE OF TRIPLE-WELL MODEL 



FIG. 1: (color online) Poincare section of the classical (mean- 
field) triple-well Hamiltonian with conjugate variables {Ji,0i) 
and {J 2 , 02 ) at 62 — 0, O 2 < 0, c = 1.25, E ss 0.708. ‘-I-’ rep¬ 
resents a state in the central regular region and represents 
a state in the chaotic sea. 



FIG. 2: (color online) Plot of |(p(a;i, * 2 , * 3 ; t)P for the quan¬ 
tum triple-well model with N = 40. Two axes are xi € [0,1] 
and X 2 € [0,1] ( 2:3 = 1 — Xi — 12 ). Red regions are of larger 
|(pp. (a) Initial state corresponding to the mean-field state 
denoted by in Fig.[H (b) the ‘-I-’ state after evolving dy¬ 
namically t = 14.5; (cjinitial state corresponding to the 
state in Fig.j^ (d) the state at t = 14.5. 


tial [331 . Under tight-binding approximation, the second- 
quantized Hamiltonian is (as a specific case of Eq. (0) 

1 c ^ 

- 2 + ajalaA , (15) 

l<ij<3 i=l 


We now illustrate our results with an example. 
Consider a BEC in a ring-shaped triple-well poten- 


where c is the on-site interaction strength. For this sys¬ 
tem M = 3. Its corresponding nonlinear mean-field EOM 
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is 


i 


_d 

dt 


( V'l 
V'2 

\ V'3 


/civ^ip -1/2 -1/2 \ /V'i\ 
-1/2 c|V'2p -1/2 U 2 
V -1/2 -1/2 dV'alV yV's/ 


Shown in Fig.is a Poincare section of the above mean- 
field dynamics, where two kinds of motion are evident: 
the central regular region is surrounded by a chaotic 
sea. The conjugate variables used in plotting Fig.[^ are 
(Ji,0i), (J2,02), which are defined as Ji = IdiP — IdaPj 
'h = IdaP, = argd 2 - argi/i, O 2 = 2 argV' 2 - argi/i- 


arg da- 

The quantum dynamics of this model can also be com¬ 
puted rather easily. The evolution of \ip{x)\‘^ in PNPS is 
plotted in Fig.[^ where two types of quantum dynamics 
are clearly observed. In Fig.[^ (a, b), an initial coher¬ 
ent state, which is a gaussian-like wavepacket in PNPS, 
shows no significant expansion or distortion during dy¬ 
namical evolution. In Fig.[^ (c, d), the situation is dras¬ 
tically different: a similar-looking initial coherent state 
expands and becomes dramatically distorted after a cer¬ 
tain time. The difference is caused by the fact that the 
initial state in Fig.[^(a) corresponds to a mean-field state 
in the regular region in Fig.while the one in Fig.[^ (c) 
corresponds to a mean-field state in the chaotic region. 

It is obvious that the mean-field theory cannot de¬ 
scribe the dramatic quantum dynamics shown in Fig.[^ 
(c, d). Such a failure or breakdown of the mean-field 
theory due to rapid decoherence has long been noticed 
in literature P5H5T] . In Ref. [50], a remedy was tried 
unsuccessfully to bridge the gap between the mean-field 
theory and the exact quantum theory. In this work we 
have shown that there exists a general time scale r/j in 
terms of Lyapunov exponent and number of bosons be¬ 
yond which the mean-field theory fails. In the following, 
we shall introduce a quantum fidelity to distinguish the 
two types of quantum dynamics shown in Fig. [^without 
using mean-field formalism, and confirm the time scale 
T/j numerically. 


A. Quantum Fidelity 

To quantify the loss of coherence in the quantum evo¬ 
lution as shown in Fig.j^ (d), we introduce the following 
quantum fidelity Tq for one-particle reduced density ma¬ 
trix (RDM) p and y: 

■^g(/5,x) = (17) 

For a quantum state |'I'(t)), its one-particle RDM can be 
explicitly written as 

I*) 01 ■ (18) 

ij 


There are three reasons to use this quantum fidelity: 


1) Experimentally we are often interested in the one- 
particle RDM. 

2) It allows us to define coherence C: 

C{p)=J^q{p,p), (19) 

where p is the one-particle RDM for Idt). The coherence 
C can quantify how coherent the state Idt) is: C{p) = I if 
and only if is a coherent state as in Eq. (§. 

3) It returns to the mean-field fidelity for coherent 

states, i.e., Fq{p,x) = = I(</I0)P if P, X are 

one-particle RDM for coherent states and 1$)^, and 

(f) are mean-field states of and |$)^ (see discussion 
under Eq. ([^). Therefore, before the Ehrenfest break¬ 
down Jjj essentially captures mean-field characteristics, 
especially the Lyapunov exponent, which distinguishes 
regular and chaotic mean-field trajectories. 


B. Numerical Results 


The numerical simulation aims at verifying our theo¬ 
retical understanding as discussed: for a coherent initial 
state, at the beginning the mean-field dynamics agrees 
with the quantum evolution, producing even the same 
growth of discrepancy between states; however, long-time 
exponential growth is not allowed by quantum mechan¬ 
ics, so there exists an Ehrenfest time Th beyond which 
the mean-field and quantum correspondence fails. Such 
a failure is due to the decoherence of quantum states; 
the breakdown time is given in Eq. 0 . In the fol¬ 
lowing we provide numerical evidences for our theortical 
understanding. 

We choose a coherent initial state |'k(t = 0))^ with one- 
particle RDM p{t = 0), whose corresponding mean-field 
state is |^(t = 0)). Then we slightly perturb the mean- 
field state into = 0 )), and generate the corresponding 
coherent state |^(t = 0))^ and RDM p{t = 0). Next we 
observe the evolution of quantum fidelity between these 
two states, which allows us to calculate the Lyapunov 
exponent. Of course, Itpit)) and |-0(t)) evolve according 
to the mean-field equations Eq. (16), |4'(t)) and |^(t)) 
evolve according to the quantum Hamiltonian in Eq. (|^ , 


p{t) and p{t) are obtained from |4'(t)) and |4'(t)), respec¬ 
tively. 1—J^q{p{t), p{t)) and I—J/nf('0(f)) V’(f)) are shown 
in Fig.[^(a), where we see that the mean-field fidelity J/if 
coincides with J/j for small t, as expected. 

However, we also observe in Fig. ^ (a) that there is 
an Ehrenfest time Th, when JFq and Jynf start to visi¬ 
bly disagree. Cases for different N and 7 are plotted 
in Fig.[^ (a), where we can see that as N increases or 7 
decreases, Th gets longer. This qualitatively agrees with 
the scaling of the Ehrenfest time. And in Fig. i (b), it 
is observed that although Th is different for different N 
and 7 , Th is approximately the time when the coherence 
C{p{t)) drops below 98%. This confirms our understand¬ 
ing that the failure of correspondence between the mean- 
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FIG. 3: (color online) (a) Quantnm and mean-field fideli¬ 
ties. Solid lines are 1 — (V’(^)) V’(^)); dashed lines are 1 — 

^q{p{i), Kt))^<i{ki),kt)) for N = 80; dot¬ 
ted lines are 1 - Fq{p{t), p{t))/Fq{p{t),p{t))Tq{p{t), p{t)) 

for N = 40. Fq is normalized to better show the correspon¬ 
dence. (b) Coherence C{p{t)). Curves show the decay of co¬ 
herence of quantum many-body states in (a). In the simula¬ 
tion, c = 1.25, E ~ 0.708, 6 i = tt, 02 = 0, p{t) and p{t) are 
the RDM of qnantum states, whose corresponding mean-field 
states are |'!/>(t)) and respectively. Wip — '<p\\t^o ~ 10 “"^. 

The lower set of lines in (a) and the corresponding upper set 
of lines in (b) are for the integrable case Ji = 0; the npper 
set in (a) and the corresponding lower set in (b) are for the 
chaotic case Ji = 0.5. 


field and quantum descriptions is the result of decoher¬ 
ence of quantum states. 


Based on such understanding, we can quantitatively 
define the Ehrenfest time in this example as the time 
when the coherence C{p{t)) drops below 98%. Examples 
of decay of C{p{t)) is illustrated in Fig.ffl (a), where the 
Ehrenfest time Th is measured when C(p{t)) drops be¬ 
low the dashed line. By varying N and c, we verify the 
relation Eq. (14), which leads to Eq. ([^, in Eig. Kb). 
A linear fitting between In A(t)/A(0) -I- Int and IniV is 
found with a constant slope A -I- ^ « 0.6 (see Eq. (|1^), 
suggesting Ci « 0.6 in Eq. ([^. Note here yt is replaced 
by In A(t)/A(0) for numerical convenience. 


FIG. 4: (color online) (a) Time evolution of coherence C{p{t)) 
for the integrable Ji = 0 and the chaotic Ji = 0.6 trajecto¬ 
ries. c = 1.25, E « 0.708, 01 = TT, 02 = 0 and time step 
in simulation is 10~®. C of integrable Ji = 0 (the solid line 
in the top) remains high while C of chaotic Ji = 0.6 (other 
solid lines) drops quickly {N = 20, 30,..., 80 from left to 
right). The dashed line is C = 98%. (b) Linear fitting of 
In A(t) /A(0) -f In t -I- In K = Ci In A for the chaotic initial state 
in (a) with c = 1, 1.25 and 2, N from 20 to 80. Data points 
are calculated when C drops to 98%. The slope ci « 0.6 is 
found independent of c or 7 . 


IV. CONCLUSION 


In sum we have answered an intriguing question — 
when does mean-field approximation of a dilute Bose gas 
remain valid as the system evolves? Our answer is the 
mean-field dynamics breaks down at the Ehrenfest time 
Th = (ci/ 7 ) In TV. The study is facilitated by introducing 
particle number phase space, where one can see easily 
that the correspondence between many-body quantum 
dynamics and mean-field dynamics is similar to the usual 
quantum-classical correspondence. 

As N can be varied in BEG experiments, it is now 
possible to experimentally measure the logarithmic be¬ 
havior of the Ehrenfest time. One can compare physical 
observables in the experiment with their theoretical 
mean-field values, and measure the Ehrenfest time 
when their discrepancy exceeds a threshold. BECs with 
unstable or chaotic mean-field descriptions are suitable 
for such experiments; for example, spinor BECs |351136j 
may be a good candidate system. 
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Appendix A: Quantum EOM in PNPS and Its 
Mean-fleld Approximation 

For the Hamiltonian in Eq. (|^, the Schrodinger equa¬ 
tion in PNPS reads {h = 1) 

idt(p{x-t) = 

i 

+ '^U^iNxi{xi - €)ip{x;t) 

i 

+ ^ Xi{xj + e)ip{x + ]t) 

i¥=3 

+ '^U^jNxiXjip{x]t), (Al) 

i¥=3 


where di = and ki{x;t) is the local wavevector of 
wavefunction ip at {x;t), as discussed in condition (ii). 
The argument of all fcj and p is (a;; t) and omitted. 

Now we assume a 5-function solution as in Eq. By 
equalling the coefficients before 5, di6 and the derivatives 
of coefficients before 6 (which is necessary to reflect the 
plane-wave phase structure) on both sides, keeping finite 
terms in the large N limit, we obtain 


dtx° = 2 Im y x°x° exp z(fc° - k°)e (A3) 

3 


where e = 1/A, e®-' is an M-dimensional vector e/' = 
—Sik + Sjk, k = 1,2,..., M. 

We are especially interested in the dynamics of a nearly 
coherent state. With conditions (i) and (ii) in Sect. II B 
and N —3 oo, Eq. (Al) becomes 


idtp = H°^Nxip -b UuNxjjx, - e)p 

i i 

-I- y^ H'I^NyJxjXjp exp[z(fcj — ki)f\ 

i¥^3 

+ d t exp[i(A:j - ki)e] 


... ^^2]jxj 


+ - di)\p\]^^exp[i{kj - ki)e] 

+ y^ UjjNxjXjp + o(l), (A2) 

i^3 


dtk°e = -Re 77° J ^ exp i{k° - k°)e 
—2 Re y^ UijX° , 


(A4) 


where the argument t of all x/ and 7/ is omitted for 
brevity. Lengthy but straightforward calculations will 
verify that Eqs. (|A3|) and (|A4|) are equivalent to Eq. ([To|, 


which is same as the mean-field EOM for the one particle 
RDM p. 
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